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We report results from a search for neutral Higgs bosons produced in association with b quarks 
using data recorded by the DO experiment at the Fermilab Tevatron Collider and corresponding to 
an integrated luminosity of 7.3 fb _1 . This production mode can be enhanced in several extensions 
of the standard model (SM) such as in its minimal supersymmetric extension (MSSM) at high tan B. 
We search for Higgs bosons decaying to tau pairs with one tau decaying to a muon and neutrinos 
and the other to hadrons. The data are found to be consistent with SM expectations, and we set 
upper limits on the cross section times branching ratio in the Higgs boson mass range from 90 to 
320 GeV/c 2 . We interpret our result in the MSSM parameter space, excluding tan B values down to 
25 for Higgs boson masses below 170 GeV/c 2 . 
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In contrast to the standard model (SM), where only 
one Higgs boson doublet breaks the SU(2) symmetry, 
there are two Higgs boson doublets in the minimal super- 
symmetric standard model (MSSM) [1] . This leads to five 
physical Higgs bosons remaining after electroweak sym- 
metry breaking; three neutrals: h, H, and A, collectively 
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denoted as <f>, and two charged, H ± . At the tree level, 
the mass spectrum of the Higgs bosons is determined by 
two parameters conventionally chosen to be tan/3, the 
ratio of the two Higgs doublet vacuum expectation val- 
ues, and Ma, the mass of the pseudoscalar Higgs boson 
A. Although tan /3 is a free parameter in the MSSM, 
large values (tan/3 > 20) are preferred. The top quark 
to bottom quark mass ratio suggests tan 8 « 35 [2] , and 
the observed density of dark matter also points towards 
high tan /3 values [3]. At high values of tan/3, two of the 
neutral Higgs bosons (A and h or H) are approximately 
degenerate in mass. They share similar couplings to 
quarks, enhanced by tan /3 compared to the SM couplings 
for down- type fermions, while the couplings to up- type 
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fcrmions are suppressed. The enhancement of couplings 
to down- type fermions has several consequences. First, 
the main decay modes of this Higgs boson pair are 4> — > bb 
and (j) — V tt with branching ratios £>(</> — >• bb) w 90% and 
B(4> —¥ tt) « 10%, respectively. Their production in 
association with 6 quarks is enhanced by approximately 
tan 2 (3 compared to the SM, which could make this pro- 
duction rate measurable at a hadron collider. 

Experiments at the CERN e+e~ Collider (LEP) ex- 
cluded MSSM Higgs boson masses below 93 GeV/c 2 [4]. 
The CDF and DO collaborations at the Tevatron 
extended the exclusion to higher masses for high 
tan j3 [5-9] . More recently, similar searches were per- 
formed at the LHC [10]. In this letter, we present a 
search for the process pp — >■ <fib — >• tt6 where one r 
lcpton (denoted t m ) decays via t — > [iv^v T and the 
other (denoted r/j) decays hadronically. This mode is 
complementary to the inclusive <p — ► TT [5, 6] and the 
<f>b — > bbb [8] searches. This is because in the former, 
the presence of b quark(s) in the final state significantly 
decreases the Z boson background, while the latter 
has a larger branching ratio but suffers from a large 
multijet background and is more sensitive to the MSSM 
parameters. This result is built on, and supersedes, 
our previous result based on 2.7 fb _1 of integrated 
luminosity [9] . In addition to the increase in luminosity, 
the sensitivity is improved by a refined treatment of 
systematic uncertainties, higher-performance signal to 
background discriminants and a higher trigger efficiency. 

The data considered in this analysis were recorded by 
the DO detector, described in [11], and correspond to 
an integrated luminosity of 7.3 fb _1 [12]. Events were 
recorded using a mixture of single high-pr muon, jet, tau, 
muon plus jet, and muon plus tau triggers. A data sample 
of Z — > T^Th is employed to measure the efficiency of this 
inclusive trigger approach with respect to single muon 
triggers. This has been validated in Z(— > T^Th )+jets 
events. The overall trigger efficiency ranges between 80% 
and 95%, depending on the kinematics and on the de- 
cay topology of the hadronically decaying r. We rely on 
all components of the DO detector: tracking, calorime- 
try, and the muon system. Muons are identified from 
track segments reconstructed in the muon system that 
are spatially matched to reconstructed tracks in the in- 
ner tracking system, and muon system scintillator hits 
must be in time with the beam crossing to veto cos- 
mic muons. Hadronic r decays are reconstructed from 
energy deposits in the calorimeter [13] using a jet cone 
algorithm with radius = 0.3 [14]. They are required to 
have associated tracks. The r candidates are then split 
in three different categories which roughly correspond to 
one-prong r decay with no ir° s (r^ type 1), one-prong de- 
cay with 7T° s (t/j type 2), and multiprong decay (r^ type 
3). In addition, we use a neural-network-based Th iden- 
tification (NN r ) to separate quark and gluon jets from 



genuine hadronic r decays [13]. The NN T is based on 
shower shape variables, isolation variables, and correla- 
tion variables between the tracking and the calorimeter 
energy measurements. We require NN T > 0.9 (0.95 for 
i~h type 3) which has an efficiency around 65% while re- 
jecting w 99% of jets. Jets are identified as clusters of 
energy in the calorimeter reconstructed with the mid- 
point cone algorithm [14] with radius = 0.5. Jet recon- 
struction and energy calibration are described in [15]. 
All jets are required to pass a set of quality criteria 
and to have at least two reconstructed tracks originat- 
ing from the pp vertex matched within Ai?(track, jet- 
axis) = (At;) 2 + (A95) 2 < 0.5 (where r\ is the pseudo- 
rapidity [16] and tp the azimuthal angle). A neural net- 
work 6-tagging algorithm [17] (NNb), with lifetime-based 
information involving the track impact parameters and 
secondary vertices as inputs, is used to identify jets from 
b quarks. The missing transverse energy, fi T , used to 
infer the presence of neutrinos, is reconstructed as the 
negative of the vector sum of the transverse energy of 
calorimeter cells with |r/| < 3.2. It is corrected for the 
energy scales of all reconstructed objects. 

The leading order (LO) event generator pythia [18] 
is used to generate <f>b production in the 5-flavor scheme, 
gb — > (f>b. To correct the cross section and the event kine- 
matics to next-to-leading order (NLO), we use MCFM [19] 
to compute correction weights as a function of the lead- 
ing b quark px and r\ in the range p T > 12 GeV/c and 
|7i fc | < 5. The dominant backgrounds to this search are 
the production of Z+jets, it and multijets (MJ), the lat- 
ter being estimated from data. We also consider VF+jets 
and diboson (WW, WZ and ZZ) production. Diboson 
events are simulated with PYTHIA while Z+jets, VF+jets, 
and tt samples are generated using ALPGEN [20] with 
pythia for showering and hadronization. tauola [21] 
is used for the decay of r leptons; b hadron decays are 
modeled with evtgen [22]. The generated samples are 
processed through a detailed simulation of the DO de- 
tector based on GEANT [23]. The output is then com- 
bined with data events recorded during random beam 
crossings to model the effects of detector noise and pile- 
up energy from multiple interactions and different beam 
crossings. Finally, the same reconstruction algorithms as 
for data are used on the simulated events. Corrections to 
the simulation are derived from data control samples and 
applied to object identification efficiencies, energy scales 
and resolutions, trigger efficiencies, and the longitudinal 
pp vertex distribution. Signal, tt, and diboson yields are 
determined from the product of the acceptance and de- 
tector efficiency (both determined from the simulation) 
multiplied by theoretical cross section times luminosity. 
For the dominant Z — > tt background, the simulation is 
corrected by comparing a large sample of Z — > /i/i events 
in data and in the simulation. This correction, measured 
in each jet multiplicity bin as a function of the <p* event 
variable [24], leading jet 77, and leading 6-tagged jet NNb, 
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D0, L = 7.3 fb' 1 data D0, L = 7.3 fb 1 -— data D0, L = 7.3 fb 1 — ~ da,a 




FIG. 1: (a) Mh at distribution in the Pretag selection, (b) T>f for a Higgs boson mass of 110 GeV/c 2 . (c) T>f for a Higgs boson 
mass of 180 GeV/c 2 . The predicted signal is shown assuming the MSSM scenario described in the caption of Table I. 



affects both the normalization and the kinematic distri- 
butions. For the VF+jets background, the muon predom- 
inantly arises from the W boson decay while the hadronic 
r candidate is faked by a jet. While this background is 
estimated from the simulation, it is normalised to data 
using a W(— > /i^)+jets control sample. 
We define a background-dominated sample, named Pre- 
tag in the following, to ensure our background modeling 
is correct. We select events with one reconstructed pp 
vertex with at least three tracks, exactly one isolated 
muon (/x), exactly one reconstructed hadronic tau (t/j), 
and at least one jet. The muon is required to have a 
transverse momentum pj, > 15 GeV/c, < 1.6, and to 
be isolated in the calorimeter and in the central tracking 
system, Ai?(/i,jet) > 0.5 relative to any reconstructed 
jet. The Th candidate must satisfy pJ^ > 10 GeV/c, 
\p Th | < 2.0, Afl(rfi, p) > 0.5 relative to any muon, and Th 
tracks must not be shared with any reconstructed muons 
in the event. We also require the distance along the beam 
axis between Th and p Az(r/ l ,/i) < 2 cm. Selected jets 
have ^ > 15 GeV/c, |?,j ct | < 2.5, AR(jet,T h ) > 0.5. 
In addition, we require Th and p, to have an opposite 
electric charge (OS) and a transverse mass Mr(p, $t) < 
60 GeV/c 2 (100 GeV/c 2 for r h type 2). The transverse 
mass of N reconstructed objects is defined as: 

M T (0 1 ,..,0 N )= I Pt -V°t ■ [1 -cosA^Oi.O,-)], 

y Oi,Oj 

where A<^(0,-, Oj) is the azimuthal angle between objects 
Oi and Oj . Most of the M J background is removed by the 
requirement 2?mj > 0.1 (0.2 for 77, type 3) where X>mj is 
a multivariate discriminant described below. Finally, to 
improve the signal to background ratio, we select a more 
restrictive 6-tagged sample by demanding at least one 
jet to have NNf, > 0.25. This 6-tag requirement has an 
efficiency of 65% for a probability of misidentifying a light 
parton jet as a b jet of 5%. Table I shows the predicted 
backgrounds, observed data yields, and expected signal 
yields in the pretag and 6-tagged samples. 

The MJ background is estimated from control data 



TABLE I: Expected background yield, observed data yield, 
and expected signal yields for the two selections described in 
the text with systematic uncertainties. The signal yields are 
given for the m™ ax scenario (p = +200 GeV and tan/3 = 40). 





Pretag 


b-tagged 


Z+jets 


2237.7 ± 123.5 


217.5 ± 16.8 


it 


225.6 ± 38.7 


182.6 ± 32.2 


MJ 


225.0 ± 39.6 


28.4 ± 4.8 


Other 


451.8 ± 18.6 


47.6 ± 3.0 


Total background 
Data 

Signal T7i0 = 110 GeV/c 2 
Signal = 180 GeV/c 2 


3139.9 ± 154.0 
3236 
107.4 
24.0 


476.0 ± 40.2 
488 
67.8 
15.0 



samples. We define a MJ-enriched control sample with 
identical requirements as in the pretag and 6-tagged sig- 
nal samples but reversing the muon isolation criteria. In 
a dedicated MJ sample obtained by requiring /i and Th to 
have the same electric charge (SS), we measure the ratio 
of the probability for a MJ-event muon to appear iso- 
lated to the probability for a MJ-event muon to be non- 
isolated: i? iso/ ^ = 7>(/i iso |MJ)/P(/^|MJ). The depen- 
dence on i] Th , p T r£ , and leading-jet pt of -R iso /-^r- is taken 
into account. This -R; so /i^ is then applied to events in the 
non-isolated-muon sample to predict the MJ background 
in the signal samples. An alternate method is used to 
estimate the systematic uncertainty. For MJ events, we 
expect the correlation between the charge of p, and Th to 
be small. Therefore, we use a data sample that has the 
same selection as the 6-tagged sample except that p, and 
Th are SS. We subtract from this MJ-dominated SS sam- 
ple the residual contribution from other SM backgrounds. 
The number of MJ events in the OS signal sample is ob- 
tained by multiplying the SS sample yield by the OS:SS 
ratio, 1.07 ± 0.01, determined in the non-isolated-muon 
sample. The difference in normalization between the two 
methods is taken as a systematic uncertainty on the MJ 
contribution. 

To further improve the signal to background discrim- 
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FIG. 2: (a) Model independent cross section times branching ratio limit as a function of (b) tan/? vs Ma limit in the 
MSSM m™ ax scenario, and (c) in the MSSM no-mixing scenario. 



ination, we use multivariate techniques. A first neural 
network X>mj is used to separate MJ background from 
the signal. Two 2?mj discriminants are trained, one for 
types 1 and 3, and another for Th type 2. They are based 
onp£, p T T h , $ T , \A<p(n,r h )\, H T = Ejets^r , M T (AUO) 
(where the sum is performed over all objects), M^ at , and 
M co i. The quantity Mh a t is defined as 

A/hat = \/ (B^ ~ P^ h + $ T ) 2 ~ |P*T h +PT + ^Tl 2 , 

where £^ Th is the energy of the r^/i system, and p^ Th 
is its momentum along the beam axis. It represents the 
minimal center-of-mass energy consistent with a di-tau 
resonance decay. The quantity M co \ is the jiTh invari- 
ant mass assuming neutrinos are emitted along the r de- 
cay axis [25]. To address the it background, we train a 
neural network T> t j to discriminate against signals built 
from samples simulated at three consecutive Higgs boson 
masses, in order to increase the signal statistics. It is 
constructed from the variables \Aip(n,Th)\, \A(p(fi, t T )\, 
H T , H T +p T T h +p T , $ T , M T {K\\0), M T (ji,$ T ), M hat , 
M co i, At = {p T —p T h )/p T h , and A/j ets , the total number of 
jets in the event. Finally, for events satisfying V t j > 0.1, 
we form a likelihood discriminant T>f which uses as input 
£> M j, £> t i, NN b , and M hat . 

Systematic uncertainties are divided in two categories: 
those affecting only the normalizations and those also 
affecting the shapes of T>f distributions. Those affecting 
the dominant -Z+jets background modeling are evaluated 
with Z — > /i/i samples: Z+jets (3.2%) and Z+6-tagged 
jets (5%) normalizations, inclusive trigger efficiency (3%) 
which also affects all other simulated processes, Z boson 
kinematics (1%) which is shape-dependent. For non-Z 
boson and non-MJ backgrounds, we consider the uncer- 
tainties affecting the normalization: luminosity (6.1%), 
muon reconstruction efficiency (2.9%), reconstruc- 
tion efficiency [(4-10)%], single muon triggers efficiency 
(1.3%), ti and diboson cross sections (11% and 7%), and 
the uncertainties affecting the shape of T>f. jet energy 
calibration (~ 10%) and 6-tagging (~ 4%). The en- 
ergy scale, and jet identification efficiencies have a negli- 



gible effect. The MJ background systematic uncertainties 
range from 10% to 40%. 

The predicted background, signal, and data distribu- 
tions of Mhat and X>/ discriminant are shown in Fig. 1. 
The T>f distributions are used as input to a signif- 
icance calculation using the modified frequentist ap- 
proach [26, 27]. We do not observe any significant ex- 
cess over the expected background. We first set model 
independent limits (assuming the Higgs boson width is 
negligible compared to the experimental resolution) at 
the 95% C.L. on the signal cross section times branching 
fraction as a function of the Higgs boson mass; these are 
shown in Fig. 3(a). These limits are then translated into 
the tan /3, Ma plane for two MSSM benchmark scenar- 
ios [28]: the m™ ax and no-mixing scenarios. The MSSM 
to SM signal ratio as well as the Higgs boson width are 
calculated with the feynhiggs program [29]. In this 
interpretation, we further include systematic uncertain- 
ties on the signal production cross section (15%) [8]. We 
also take into account the Higgs boson width using the 
method described in [8]. Figures 3(b) and (c) present the 
limits for the two scenarios with the higgsino mass pa- 
rameter fi = +200 GeV. Numerical results and limits in 
other MSSM scenario are presented in [30]. We exclude 
a substantial region of the MSSM parameter space, es- 
pecially at low Ma, and set the most stringent limit to 
date at a hadron collider, when using this final state. 
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FIG. 3: Limits on tan /3 vs Ma for different benchmark scenarios: (a) m™ ax with 
M = -200 GeV. 



-200 GeV, (b) no-mixing with 



TABLE II: Expected and observed upper limits on tan /3 as a function of Ma in four MSSM benchmark scenario. 



M A 




m,h max 




mi max 




no mixing 




no mixing 


(GeV/c 2 ) 


M 


= -200 GeV 


A* 


= +200 GeV 




= -200 GeV 




= +200 GeV 




Obs. 


Exp. 


Obs. 


Exp. 


Obs. 


Exp. 


Obs. 


Exp. 


90 


28.9 


21.6 


28.8 


21.3 


29.0 


21.6 


28.8 


21.2 


100 


32.0 


23.5 


32.0 


23.5 


31.9 


23.3 


31.7 


23.2 


110 


33.8 


24.7 


33.6 


24.4 


33.6 


24.9 


33.5 


24.4 


120 


24.8 


26.6 


25.0 


26.7 


25.0 


27.0 


24.9 


26.8 


130 


26.7 


27.8 


26.3 


27.5 


26.7 


27.9 


26.6 


27.6 


140 
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29.9 


28.1 


29.3 
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29.5 
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29.4 
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25.3 


31.5 


25.3 
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25.6 


31.6 


25.3 


31.3 


160 


27.5 


32.8 


27.3 


32.6 


27.2 


32.7 


27.4 


32.4 


170 


30.1 


36.3 


29.8 


36.0 


30.0 


36.2 


29.5 


35.8 


180 


39.0 


41.5 


38.4 


40.8 


38.4 


40.9 


38.4 


40.7 


190 


41.9 


44.1 


41.6 


43.4 


41.6 


43.8 


41.5 


43.6 


200 


46.9 


48.6 


45.8 


47.1 


46.6 


48.2 


45.5 


47.4 


210 


52.8 


53.9 


51.8 


53.3 


52.3 


53.5 


51.9 


53.1 


220 


57.2 


58.9 


56.7 


58.0 


57.0 


58.5 


56.5 


57.8 


230 


61.8 


64.2 


60.9 


62.9 


61.9 


64.4 


61.4 


63.7 


240 


69.8 


68.7 


68.6 


67.7 


69.5 


68.6 


68.7 


67.8 


250 


75.3 


74.6 


75.5 


74.8 


74.7 


74.1 


72.9 


72.0 


260 


79.2 


78.8 


81.0 


80.4 


80.7 


80.3 


79.3 


78.8 


270 


91.1 


86.4 


93.5 


88.8 


92.5 


87.4 


91.1 


86.5 


280 


96.1 


91.4 


99.2 


94.1 


97.2 


92.9 


96.8 


92.0 



